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Abstract 

In supersymmetric scenarios with a light stop particle ti and a small mass difference to the light- 
est supcrsymmctric particle (LSP) assumed to be the lightest neutralino, the flavour changing 
neutral current decay ti — > cxi can be the dominant decay channel and can exceed the four- 
body stop decay for certain parameter values. In the framework of Minimal Flavour Violation 
(MFV) this decay is CKM-supprcssed, thus inducing long stop lifetimes. Stop decay length 
measurements at the LHC can then be exploited to test models with minimal flavour breaking 
through Standard Model Yukawa couplings. The decay width has been given some time ago by 
an approximate formula, which takes into account the leading logarithms of the MFV scale. In 
this paper we calculate the exact one-loop decay width in the framework of MFV. The compari- 
son with the approximate result exhibits deviations of the order of 10% for large MFV scales due 
to the neglected non-logarithmic terms in the approximate decay formula. The difference in the 
branching ratios is negligible. The large logarithms have to be resummed. The resummation is 
performed by the solution of the renormalization group equations. The comparison of the exact 
one-loop result and the tree level flavour changing neutral current decay, which incorporates the 
resummed logarithms, demonstrates that the resummation effects are important and should be 
taken into account. 



1 Introduction 



The Standard Model (SM) provides a very successful effective theory of particle interactions which 
is in excellent agreement with electroweak precision data. Furthermore, remarkable consistency 
and precision tests have been made in the sector of quark flavour violation. These tests and limits 
on flavour changing neutral currents (FCNC) from K, D and B studies put strong constraints on 
possible New Physics beyond the SM [1] . They forbid a generic flavour structure of New Physics 
at the TeV scale and raise the question why contributions from New Physics at ~ 1 TeV are 
strongly suppressed. A solution to this New Physics Flavour Puzzle is provided by the framework 
of Minimal Flavour Violation (MFV) [2]-[5]. It requires all sources of flavour and CP-violation to be 
given by the SM structure of the Yukawa couplings. Flavour mixing in models of New Physics is 
then always proportional to the off-diagonal elements of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix 0. In particular, the mixing of the third generation squarks with the first and second 
generation squarks is highly suppressed by small CKM quark mixing angles. 

Since the flavour structure of New Physics at the TeV scale must be non-generic, flavour mea- 
surements provide a good probe of New Physics. One of the best studied examples is given by 
supersymmetry (SUSY). In the context of flavour violation the phenomenology of the light stop 
partner ti is especially interesting. In most SUSY models a light stop quark mass arises naturally. 
Due to the large top Yukawa coupling the mixing between the weak eigenstates ti and Ir leads to a 
large mass splitting between the stop mass eigenstates. Furthermore, the large top Yukawa coupling 
in general drives the stop mass via the renormalization group equation (RGE) running to smaller 
masses, even if all squarks have a common mass at the SUSY breaking scale. In scenarios with 
a light stop which predominantly decays into a charm quark and the lightest neutralino, assumed 
to be the lightest supersymmetric particle (LSP), ii — )■ ex?, squark flavour violation can be tested 
by exploiting stop decay length measurements [7]. It has been shown, that in these scenarios light 
stops can be discovered at the LHC [8]. Finally, a light stop is also favoured by baryogenesis which 
requires a stop mass with about the top mass value or less for successful electroweak baryogenesis 
within the context of the Minimal Supersymmetric Extension of the Standard Model (MSSM) [9]. 

Minimal Flavour Violation naturally arises in supergravity models which provide flavour-inde- 
pendent scalar mass terms at a high scale like the Planck scale Mp. Some time ago the authors of 
Ref. [To] have provided an approximate formula for the flavour changing neutral current (FCNC) 
decay ii — ?• ex? by starting from a vanishing tree level ti — e — x? coupling at the Planck scale. 
The decay then proceeds via one-loop diagrams. The non-vanishing divergences, which are due to 
scalar self-mass diagrams, have been subtracted by a soft counterterm at the Planck scale so that a 
large logarithm ln(Mp/M^) remains at the weak scale chosen to be M^. The authors argue that 
in view of this large logarithm, the remaining non-logarithmic part of the one-loop diagrams can 
hence safely be neglected so that their result for the decay width takes a rather simple form. 

In this work, we perform the complete one-loop calculation of the ti — >• ex? decay in the frame- 
work of MFV. We perform the full renormalization program and keep the flnite non-logarithmic 
terms arising from the loop integrals. This allows us to study the importance of the neglected 
non-logarithmic pieces in the formula given by Ref. |10j . 

To get a reliable result, the appearing large logarithms should be resummed. In the renor- 
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malization group approach this corresponds to solving the renormahzation group equations for the 
scalar soft SUSY breaking squark masses. As has been pointed out in [1], the hypothesis of MFV is 
not renormahzation group invariant. Flavour off-diagonal squark mass terms are hence induced by 
the Yukawa couplings, so that the squark and quark mass matrices cannot be diagonalized simul- 
taneously any more and the stop state receives some admixture from the charm squark, inducing a 
FCNC between stop, charm and LSP neutralino, ii—c—Xi- From this point of view, the logarithmic 
piece of our one-loop result is equivalent to the first order in the expansion of the RGE solution for 
the squark-quark-neutralino coupling in powers of a, whereas the tree level decay calculated with 
the FCNC coupling includes the resummation of the large logarithms. The comparison of the two 
decay widths provides an estimate of the importance of the resummation of the large logarithms. 

The outline of our paper is as follows: In section [2] we present the diagrams contributing to 
the one-loop decay. We set up our notation for the squark and quark sector in section [3l The 
counterterms and the renormahzation are discussed in section [H Section [5] contains the numerical 
analysis. We conclude in section [6l In the Appendix we list our Feynman rules, the various 
amplitudes contributing to the decay and we derive the FCNC counterterm. 

2 One-loop decay 

We work in the framework of the MSSM with MFV so that all flavour changing effects with quarks 
and squarks are controlled by the quark Yukawa couplings and CKM mixing angles [3] . The decay 
of the lightest stop ti into the lightest neutralino Xi aiid a charm quark c. 



is then mediated at the one-loop level. We consider scenarios where the light stop ii is the next- 
to-lightest super symmetric particle (NLSP) and the lightest neutralino is the LSP. The process is 
built up by the stop and charm self-energies and the vertex diagrams, cf. Fig. [TJ Note, that in our 
calculation we set 



Therefore in the ti self-energies we have only non-vanishing contributions for transitions into the 
left-handed charm squark cl. Those into right-handed scharm, cr, are zero for rric = 0. All dia- 
grams are mediated by charged current loops. The various diagrams which contribute are depicted 
in Fig. [2j The self-energies and vertex corrections are divergent and have to be renormalized. The 
counterterms for the squark and quark self-energies and for the vertex renormahzation are shown 
in Fig. O The FCNC vertex does not arise at tree level. Its occurrence as counterterm at one-loop 




(1) 



nic = . 



(2) 




Figure 1: Generic diagrams contributing to the loop-decay ti — > cxi- 
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level is due to the fact that MFV is not RGE- invariant, since the weak interactions affect the squark 
and quark mass matrices differently [TT] . Their simultaneous diagonalization cannot be maintained 
at higher orders so that it can consistently be imposed at a single scale only, called hmfv in the 
following. 

For the calculation of the stop decay process we define an effective interaction vertex 

T = gucik2) {FlVl + Fi^Vr) v^o{ki) , (3) 

where Uc,v^o denote the charm and neutralino spinors and ki,k2 are the four- momenta of the 
outgoing neutralino and charm quark. Fl and Fr are form factors associated with the chirality 
projectors Vl and Vr, respectively. They receive contributions F"" from the vertex diagrams, 
ptc^ptiCL from the quark and squark self-energies, and 6F'" ,6F^'^,6F^^'^^ from the vertex and the 
quark and squark wave function renormalization counterterms, 

Fl,r = [F*i^^ + 5f'^''^^ + F*" + 5F*^ + + SF'']l,r . (4) 

They are specified in section H] and the Appendices Bl-3 and C. 



3 The quark and squark sector 

For the choice of our notation the quark and squark sectors are discussed here in more detail. 
The definition of the couplings is deferred to Appendix A. We define the 3x3 unitary matrices 
U^L,R^Ud.L,R g^g ^]^Q matrices which rotate the left- and right-handed up- and down-type quark 
current eigenstates UL,R,dL^R to their corresponding mass eigenstates, u^R,d^R-, 

ut = U''^UL, u'^ = U^^UR, df = U''^dL, d'R = U''^dR. (5) 

The CKM matrix V is given by 

V = W^U'^^'^. (6) 
For the squark interaction eigenstate we define a six component vector 



q'R 



(7) 



where q'l^q'R are each a three component column vector in generation space. The squared squark 
mass matrix can be written as a 2 x 2 Hermitian matrix of 3 x 3 blocks 

It is diagonalized by a 6 x 6 unitary matrix which rotates the squark interaction eigenstates to 
their mass eigenstates q^, 

= Wq' . (9) 
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The six component column vector is defined to be ordered in mass, with qi being the lightest 
squark. Equation ([9]) can then be rewritten as 

C = Ws^qiL + Ws^+3qiR (s = 1, .., 6, i = 1, 2, 3) 

= mq'L + WRq'j,),, (10) 

where i denotes the generation index0 The rotation of the squarks by the same unitary matrices 
U''^-^ as the quarks defines the super-CKM basis. In models with non-minimal flavour violation 
the squark mass matrix is flavour-mixed in this basis, in contrast to the quark mass matrix. In 
models with MFV at the scale //mfVj however, the squarks can be rotated by U^^'^ to their flavour 
eigenstates in parallel to the quarks, and the super-CKM basis is at the same time the flavour 
eigenstate basis. Hence, suppressing generation indices, 

QL = U^^q'L, qR = U'l^q'^ . (11) 

The squared mass matrix in the flavour eigenstate basis {qL,QR)'^ then reads 

^2 ( iMl+ml)ls rnqiAq-,r)lA ^^^^ 



" » m,(A,-;ur,)l3 (M4+7n2)l 
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where = 1 /r„ = tan f3 for down- and up-type squarks. With tan (3 we denote the ratio of 
the vacuum expectation values of the two complex Higgs doublets. They are introduced in order 
to generate masses of up- and down-type fermions [12]. The parameter Aq denotes the trilinear 
coupling of the soft SUSY breaking part of the Lagrangian, fi the Higgsino mass parameter and 
niq the quark partner mass. I3 is a 3 x 3 unit matrix in generation space. The parameters Mg^ ^ 
are given by the left- and right-handed scalar soft SUSY breaking masses Mg^ ^ and the i^-terms, 

M? — M? -I- D- 
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Dq^ = Micos2(3{I;l -Qqsm^ew) 

'q^ = M| COS 2/3 Qg! 



Dq^ = M| COS 2/3 Qg sin2 , (13) 



where Ig denotes the third component of the weak isospin, Qg the electric charge, Mz the Z boson 
mass and 9]y the Weinberg angle. The squared mass matrix can be diagonalized by a 6 x 6 
unitary matrix which rotates the flavour eigenstates to their mass eigenstates, 

qT = Wstl ] =WsiqLi + Ws^+3qRi = {WLqL + WRqR)s (14) 



t 



(s,t = l,..,6, i = 1,2,3) 



Comparison with Eqs. (jlOp and (]lip shows that we can factorize the 6x3 matrices Wl^r into the 
6x3 flavour-diagonal matrices Wl^r and the 3x3 quark rotation matrices defined above, 

Wl = WlU''^ and Wr = WrW"" , (15) 

with q = u,d. The matrix W can be expressed in terms of mixing angles by 

(Wlh = iWR)^+3^ = COsOq^ , {WrU = -iWL)^+3^ = Sm9q^ . (16) 



^ We have thus decomposed the mass eigenstate squark field into left- and right-chiral interaction eigenstate squarks. 
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For the three quark generations i the relation between the flavour eigenstates qn, QiR and the squark 
mass eigenstates g"^ = (^^,^3) hence reads 



QiL COS 9q, + qm sin 9q. 
-qn sin Oq^ +qiR COS 0q^ 



(17) 



For better legibility, we suppress the generation indices from now on wherever possible, and the 
lighter and heavier squark mass eigenstates are generically called qi and q2- The mixing angles are 
then given by 

Ml -Ml 



sin 29„ 



2mq{Aq - firq) 
M? -M? 



cos 29, 



IL 



QR 



" M? -M? ' 

91 92 



and the masses of the squark mass eigenstates read 



M, 



91,2 



2 1 

g -r 2 



^4 + ^4 T - ^4)' + - /^^9)' 



(18) 



(19) 



Since the mixing angles are proportional to the masses of the quarks, the mixing is important in 
the stop sector and can drive the lightest stop mass even lighter than the top quark mass. 

4 Counterterms and Renormalization 



The quark and squark self-energies and the vertex diagrams are ultraviolet (UV) divergent and 
need to be renormalized. The quark and squark wave functions are renormalized on-shell. With 
the bare (s)quark fields q^^^ (q^^^) related to the renormalized (s)quark fields q (q) by 



1 + -6Z'J ] q 



and gg)^ = (1 + ISZ"^'""^ qL,R 



(20) 



this leads for the squarks to the following off-diagonal elements of the wave function renormalization 
constants in terms of the real part of the squark self-energy S, 

2 



2 ^ ReSst(m|) s,t = l, .., 6, s / t 



(21) 



Defining the following structure for the quark self-energy, 

we have for the corrections to the off-diagonal chiral components of the quark wave functions 
2 



(22) 



ml - 



m„ 



Re^ff{ml^) + rug^ ReS --^ (771^^. ) + nig^ ReS-J(mg^,) + mg^nig^ ReS--(m, 



4 



^Rs 



^2 



m 



m„.mn 



ReS^.^(m2^.) + mg^nig^ ReS.f (m^^) + m.^m^^. ReSf,(m^J + 



^Rs 



2 



m^.ReSg(m^,) 



iJ = 1,2,3, iy^j. 



(23) 



As for our scenario we chose the ti to be the NLSP, all our self-energies are real, and after the 
on-shell renormalization of the quark and squark wave functions we are only left with the one- 
loop vertex diagrams and the FCNC vertex counterterm. From now on "Re" will be dropped. We 
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regularize the divergences by dimensional regularization in n = 4 — 2e dimensions so that ultraviolet 
singularities appear as poles in e. We have explicitly verified that the same result is obtained with 
dimensional reduction. By exploiting the unitarity relations of the CKM matrix as well as those of 
the chargino mixing matrices U and V , defined in Appendix A, we find for the vertex contribution 
to the form factors 



an 



1 



+ ln 



A* 



+ finite terms 



loop 







(24) 
(25) 



where denotes the 't Hooft mass of dimensional regularization. We have used the short-hand 
notation 



1 



IGvr- 



r5'V2 



Zii, , Zi2 



Vcb V;. ml COS Ot 



2A/f2 r2 



(26) 



where Zii,Zi2 denote matrix elements of the Ax A Z matrix, which diagonalizes the neutralino 
mass matrix, cf. Appendix A. The further 'finite terms' in Eq. (j24p . which do not depend on In/x^, 
can be extracted from the full one-loop result for the vertex which is given in Appendix B3. We 
have introduced a generic mass for the loop particles, mioop. Note, that in the numerical analysis 
we will use the exact results with the different loop particle masses. Here, for reasons of legibility 
and also to make later contact with the result derived in Ref. [TU] we adopt the generic notation. 
The left-handed form factor is zero due to our choice of vanishing c-quark mass. 

The FCNC counterterm arises from the flavour non-diagonal part of the wave function renor- 
malization, from the renormalization of the quark and squark mixing matrices jl3H16] and the 
renormalization of the quark massesEI The renormalized squark mixing matrix is related to 
the bare by 



^in^ = {5 St + 6Wst)Wi 



tu ' 



s,t,u = 1, .., 6. 



(27) 



And similarly, the renormalized quark mixing matrices [/"^.^ *" are related to the bare quark matrices 

UUL,R{0) 



UL,B. (0) 



ik 



ij,k = 1,2,3. 



(28) 



The indices s,t,u denote the six squark mass eigenstates, and i,j,k are generation indices. We 
impose the MFV condition on the renormalized mixing matrices , and hence demand them 
to be flavour diagonal. This leads to flavour non-diagonal counterterms 6vj,5u^^-^ . Furthermore, 
as the bare and renormalized mixing matrices are unitary the counterterms must be antihermitian. 
The UV divergent part of each counterterm is determined such that it cancels the divergent part 
of the antihermitian part of the corresponding wave function renormalization matrix |14H16| , 



6w 



^{5Z^-6Z^^) 



L.R 



(29) 
(30) 



In our renormalization procedure and the definition of the counterterm we follow the same approach as in Ref. [17j . 
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The general form of the FCNC vertex counterterm depicted in Fig. [3] has been derived in Appendix 
C. For our process we have the following form factor contributions to the vertex counterterm, 



-ig ell 



]-5Z^.^ cos Ot + \5Z^, J + ^uil- cos Qt + 



0, 



(31) 
(32) 



with 



and 



V2 



Z 



(33) 



5Z, 



ct 



6Z'! r 

cz,ti 



2S,^/(m^ = 0) = 



2T!i r (m? ) 



ml 



2 

ti 



(34) 
(35) 



The squark wave function has been renormalized at p"^ = m-^ . The finite parts of the countert' 



J— — ■— — K 

erms Eqs. (|29|30p depend on the renormalization conditions. Absorbing also the finite part of the 
antihermitian wave function renormalization leads to a gauge-dependent on-shell renormalization 
scheme |15 p i6 1 [T8]. Performing minimal subtraction on the other hand is gauge-independent |19j 
and imposes the MFV condition on the MS parameters at the scale /^mfv- Ii^ the following we 
will adopt this scheme. Consequently, the result will depend on the MFV scale //mfv- The squark 
mixing matrix counterterm then read^ 



5w 



t _ 



1 Tfi - (7n? ) + S 



m 



(36) 



MS 



A gauge invariant prescription for the quark mixing matrix counterterm is given by [15| . 



[j^tiO) + 2S,^/(0)]5^ = --si;(o)j 



For the quark mixing matrix contribution to the vertex counterterm we then find 



(37) 



^e^icos^t _ 

MS 



g- 



- + In — 



For the contribution from the squark mixing we have 



with 



geii 



6w 



CLhJuS 



-ra- 



ni 



^ = g- 



cl 



2X 



2 2 

TOj — m- 



- + In^ 

£ /^MFV 



(38) 



(39) 



A 



+ + + 4(M^(t^ - 1) + M\tl) + ratA^ tan ( 



(40) 



This counterterm definition can lead to large contributions to the matrix element if mi^ ~ mg^^ . For a discussion, 
see e.g. Refs. [20] . In the scenarios of our numerical analysis we have m^^ <g f^Zi^. 
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It depends on the Higgsino parameter ^u, the soft SUSY breaking mass parameter Mj^^ including 
D term contributions, the triUnear couphng Ab, the mixing angle (3, the W boson mass and the 
pseudoscalar Higgs mass Ma- 

The contribution from 6Zl r is given by 




- +ln- 



'loop 



+ finite terms 



(41) 



As before, we have introduced a generic loop particle mass, and 'finite terms' denote further terms 
which do not depend on In p?. Their specific form can be extracted from the explicit formulae of the 
stop self-energies given in Appendix Bl. Inserting Eqs. (|34|38|39|4ip in Eq. ()3ip . the right-chiral 
part of the FCNC counterterm is then given by 

"l ml +A u2 ml +A 

e 



g6F]^ = igT 



2 

ti 



In 



+ 



■In. 



+ finite terms 



m^ 



mt 



(42) 



-ti "^CL ■■''loop 

Adding Eqs. and ([12]) and replacing F and A by Eq. (p6]) and Eq. ([l0|). respectively, we arrive 
at the following final result for the form factors, which contribute to Eq. ([3]), 

VcbVtb^b cos 6*4 



gFR 



167r2 



^V2 



Zll Zi2 



1 



2Af, 



2 „2 



m% 



[mi^ -H^ + Ai + + cjiM^itj - 1) + Mitj) + mtAb tan 



In 



[ /fMFV ] 
V<op/ 



+finite terms 
. 



(43) 



(44) 



Finally, the stop decay width in terms of the form factors Eqs. ()43l44p is given by 

2 



r(ii 



CX?) 



IGtt 



1 




(45) 



As can be inferred from F^, depending on the scale of MFV, the logarithm can become very large, 
and the decay can become important in certain regions of the parameter space, especially for large 
values of tan/3. The finite terms, which do not depend on In/z^py, are then only subleading. If 
we drop the finite terms in Eq. (|43p. the approximate result given by Hikasa and Kobayashi in 
Ref. [To] should be reproduced. In fact, for ttt-c = we can rewrite 

fi^ + 4 {M^ (4-1) + Mlt}) = Mfj^ + Ml + \a4 sin2 Ow cos 2/3 , (46) 



m 



where Mjj^ denotes the mass parameter of the Higgs doublet which couples to down-type 
fermions. With this relation the form factor Eq. (j43|) leads to the approximate result of 
Ref. [TU], if we set the MFV scale equal to the Planck scale, fJ-uFV = Mp, choose Mw as generic 
loop particle mass and neglect the finite terms. 



-^11, , Zi2 



[Mfj+Ml+Al + Ml 



I Vcb Vfb mj cos 
mtAb tan 9t] In 




(47) 
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In our full one- loop calculation of the decay width in terms of Fr, Eq. ()43p . the finite terms 
are included, and the relevance of these contributions can be checked by comparing with the 
approximate result F^/^ for the decay width in terms of the form factor F^^^ , Eq. (j47p . This will 
be discussed in section [5l 

To get a reliable result, the large logarithms of the MFV scale in the decay formula should be 
resummed. The logarithm is related to the running of the FCNC coupling of the neutralino to a 
quark and squark of different generations. We have required this coupling to vanish at the scale 
A^MFV- Minimal Flavour Violation is not RGE-invariant, however. Even though MFV is imposed 
at hmfv, at any other scale /x 7^ /^mfv a FCNC coupling will be generated through renormalization 
group evolution. The solution of the one-loop RGE for the quark and squark mixing matrices 
provides the resummation of the large In/ul^py. The coefficient of this logarithm in Eq. (I43p is 
then given by the first order in the expansion of the RGE solution for the squark-quark-neutralino 
coupling in powers of a. In the following we will call the right-handed form factor including the 
resummation effects F^^ . It is given by the FCNC coupling obtained through renormalization 
group evolution including flavour violation of the squark and quark mixing matrices, from some 
high scale down to the scale relevant for the decay process. 

5 Numerical Analysis 

The scenarios for the numerical analysis have been chosen such that they lead to a NLSP stop ti 
and a Xi LSP. The latter represents a promising dark matter (DM) candidate in the MSSM |21j . 
The mass difference is chosen to be small enough so that the loop mediated flavour changing decay 
^1 ~^ is dominating and can compete with the four-body decayj^ into the LSP, a 6-quark and 
a fermion pair |23] . 

ii ^ Xibff . (48) 

Such scenarios can be consistent with electroweak baryogenesis [9] and also with Dark Matter 
constraints [24]. The mass spectra and mixing angles have been calculated with the spectrum 
calculator SPheno [25] and compared to SOFTSUSY [26]. Both codes include the option to perform 
two-loop RGE running with and without the inclusion of flavour violation and both support the 
SUSY Les Houches accord |27j . Within this accord the gauge and Yukawa couplings as well as 
the soft SUSY breaking mass parameters and trilinear couplings are given out as DR running 
parameters at a scale Q, which we have chosen to be the scale of electroweak symmetry breaking 
(EWSB). We have verified that the calculation of the decay width leads to the same result if 
we apply dimensional reduction instead of dimensional regularization, so that the DR running 
parameters can be used. The mixing matrix elements and the SUSY particle pole masses have 
been taken at the scale of EWSB as well. The SM parameters have been chosen as Mz = 91.187 
GeV, a-£^{Mz) = 127.934, af^{Mz) = 0.1184, mp(mfe) = 4.25 GeV, Mf°^° = 173.3 GeV and 
^poie ^ ^ j^j j^^^g chosen the CKM matrix elements as \Vtb\ = 0.9993 and \Vcb\ = 0.04. 

In order to ensure MFV, for all three generations a common mass parameter Mg^ for the soft SUSY 
breaking masses of the SU{2) doublet has to be introduced at the scale /i = ^mfVi so that the 

^Scenarios where 2-body decays at tree level are forbidden for the lightest stop quark and where the loop induced 
flavour changing decay competes with 3-body decays have been discussed in [22] . 
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up- and down-type squark mass matrices can be simultaneously flavour-diagonal. We work in the 
framework of a MFV MSSM defined at the GUT scale in terms of a small number of parameters. 
They are given by common soft SUSY breaking scalar and gaugino mass terms, Mq and M1/2, 
a common SUSY breaking trilinear coupling Aq, the ratio of the two vacuum expectation values 
tan (3 and the sign of the Higgsino parameter /i. 



5.1 Analysis for /imfv ~ 10^^ GeV 

We first investigate two mSUGRA scenarios with soft-breaking terms at the GUT scale Mqut ~ 
10^^ GeV, which is identified with the MFV scale. All soft SUSY breaking parameters are family 
universal. The boundary conditions at /Umfv = Mgut are 

(1) Mo = 200 GeV M1/2 = 230 GeV Aq = -920 GeV 
tan/3 = 10 sign(^) = + 

(2) Mo = 200 GeV M1/2 = 230 GeV Aq = -895 GeV ^ ' 
tan /3 = 10 sign(/i) = -|- . 

The second scenario has a larger ii — Xi mass difference compared to scenario (1), whereas the 
mass difference between ii and the lightest chargino xt smaller. Since the 4-body decays are 
dominated by the chargino exchange diagram |23] . in scenario (2) the 4-body decays should be more 
important leading to a smaller branching ratio of the flavour changing decay. The GUT scale is 
given by ~ 2.3- 10^^ GeV. The masses are obtained by RGE evolution from the GUT scale down to 
the electroweak scale. The running is performed at two-loop order without the inclusion of explicit 
flavour violation in the squark sector. The obtained masses are 

(1) mr = 104 GeV m^o = 92 GeV m-+ = 175 GeV 

(2) mr = 130 GeV = 92 GeV m-+ = 175 GeV . ^ ^ 

For these scenarios the partial stop decay width into charm and neutralino, calculated with the full 
one-loop formula, is compared to the approximate result. For the latter, we take Mw as generic 
loop particle mass, cf. Eq. (I47p . The widths and form factors are given in Table 1. They have been 
obtained with the program SUSY-HIT [28], where the full one-loop formula for the flavour changing 
stop decay has been implemented. As can be inferred from the table, the exact and approximate 







1 T-|l-l00p| 

\ 


rfi/K[GeV] 


iFri 


Scenario(l) 


9.322 • 10^1° 


1.486 • 10~^ 


1.004 • 10-^ 


1.542 • 10"^ 


Scenario (2) 


5.862 • 10~9 


1.460 • lO^'' 


6.446 • 10"^ 


1.531 • 10"^ 



Table 1: The partial widths and form factors for the decay — )■ cXi in two MFV scenarios, calculated with 
the exact 1-loop formula, r^"'°°P, and with the approximate formula of Ref. [TU], F^/^, F^^^. 



decay width differ by 0(10)%. In fact, the finite terms extracted from the one-loop formula turn 
out to contribute with ~ 3 — 5% to Fr, Eq. (|43l) . This difference leads to the 10% effect in the 
decay width. The difference in the branching ratio BR{ii — t- cxi) calculated in the two approaches 
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is negligible, however. We note, that in the first scenario the partial width is ~ 6 times smaller 
than in the second scenario due the smaller ii — x? mass difference and hence reduced phase space. 

For the calculation of the branching ratios, also the partial width for the ti decay into u-quark 
and neutralino, ii — t- nx?, as well as the 4-body decay width are needed. The former is suppressed 
by 2 orders of magnitude compared to the cxi final state due to the small CKM matrix element 
\Vub\ ~ 0.003 which enters quadratically in the decay width. The branching ratios are listed in 
Table 2. As anticipated, the stop 4-body decay is more important in scenario (2) leading to a 
change of the branching ratio of interest, BR(ti — t- Xi^), at the few per-cent level. 



branching ratio 


BR(ti ^ xic) 


BR(ti ^ x» 


BR(ti ^ xlbff) 


Scenario(l) 


0.9944 


0.0056 


4.587 • 10^5 


Scenario (2) 


0.9443 


0.0053 


0.0504 



Table 2: The <i branching ratios for different final states for scenario (1) and (2). 



As stated before, the large logarithms in the decay formula should be resummed. To get 
an estimate of the importance of the resummation effects, the one-loop decay calculated in the 
framework of MFV is compared to the tree level stop decay into charm and neutralino with flavour 
off-diagonal elements in the squark mixing matrix. They are the result of the small flavour off- 
diagonal entries introduced in the soft-breaking terms through RG evolution including the complete 
flavour structure of the different flavour matrices, from the scale of MFV down to the scale of EWSB. 
The input parameters for the decay formula are taken from SPheno [250. In the flavour violating 
case, denoted by FV in the following, where no flavour-eigenstates exist any more, the lightest 
up-type squark state ui has been identified to correspond to ii. The scenarios in Eq. (j39|) have 
been chosen such that m^^ ~ m^^. The x? a-nd Xi^ masses are almost unchanged. The form factor 
of the tree level decay is given by the right-handed part of the FCNC ui — c — x? coupling, 

F^ = -^V2^^tw + ^)m),,^, (51) 

with the squark mixmg matrix Wl defined in Eq. (fTO|) . This leads to the partial decay width 

The form factors and partial widths are shown in Table [3l As can be inferred from the table, 
there is a factor ~ 4.4 between the right-handed form factor calculated at the one-loop level in the 
MFV framework and the one derived from RG evolution including flavour violation. As expected, 
resummation effects turn out to be important for a large scale ;Umfv = -^gutIiI- The partial widths, 

■'Thanks to Werner Porod who provided us with the newest SPheno version 3.0.beta56. 
^The left-handed part is negligibly small for rric = 0. 

^This result is in agreement with the discussion in Ref. where resummation effects in the coupling fi — c — Xi 
have been found to be large. 
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1 pl-loop| 




pl-loop [QgY] 


r^v [GeV] 


Scenario (1) 


1.486 • 10"^ 


3.361 • 10-5 


9.322 • 10^1° 


4.766 • 10-^^ 


Scenario (2) 


1.460 • 10-"^ 


3.306 • 10-5 


5.862 • 10-9 


3.006 • 10-1° 



Tabic 3: The right-handed form factors and partial decay widths of the lightest up-type squark into charm 
and neutralino for the MFV scenario (1-loop) and the FCNC tree level decay (FV). 

which depend quadratically on the right-handed form factor, differ by a factor ~ 20. 

For comparison we have performed the calculation with the decay spectra and mixing angles 
evaluated by SOFTSUSY. The squark mixing matrix elements agree within 10-^ accuracy with 
the results of SPheno. The mixing matrix element (Wl)mic, which enters in the form factor 
Eq. ()5ip . is 0(10-^) and differs in the two spectrum calculators. The two codes implement the 
one-loop corrections to the squark mass matrices differently. SOFTSUSY corrects only the flavour- 
diagonal entries of the squark mass matrices, while SPheno implements a full one-loop calculation, 
so that differences in the flavour off-diagonal entries are to be expected. For the SOFTSUSY 
parameter values, this results in a ratio between loop decay and tree level decay of ~ 2.7 for the 
two scenarios, compared to the ratio ~ 4.4 found with the SPheno parameter values. All in all, the 
results with both spectrum calculators show the importance of the resummation effects. 

Of phenomenological interest are the consequences of these resummation effects on the ti 
branching ratio into the charm plus neutralino final state. To quantify this, the competing 4- 
body stop decay width, calculated in the FV scenario and including tree level FCNC couplings, is 
needed. The additional FCNC contributions are expected to be small, however, due to the suppres- 
sion by CKM matrix elements. As the calculation is not available at present we defer a detailed 
comparison to a future publication. 

5.2 Analysis for /^mfv < ^gut 

In the previous section the importance of resummation effects has been discussed. With decreasing 
/xmfv and hence smaller In fJ-y^-py on the other hand the non-resummed one-loop MFV result should 
approach the resummed flavour- violating tree level result. Furthermore, we expect the approximate 
formula of Ref. [10], which is a good approximation of the exact one-loop MFV result for large 
scales, to be less good with decreasing MFV scale. In order to verify this behaviour we have chosen 
scenarios with different //mfv varied between 10^ GeV < ^mfv ^ 10-"^^ GeV. We choose the soft 
SUSY breaking input parameterj§ in each scenario such that the masses for ii and Xi remain 
almost unchanged. Consequently, the differences in the partial decay widths will not be due to 
phase space effects. Furthermore, the scenarios are constrained by the requirement that ti is the 
NLSP and x? is the LSP, so that the FC decay ii — )• ex? dominates. The relevant particle masses 
for the different scenarios vary as 

m^-^ = 105 ... 1 16 GeV and m^o = 92 ... 104 GeV , (53) 
*In our scenarios tan /3 varies between 10 and 20. 
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Q I I I I I I I I 

10^ 10^ 10** 10** 10"' IQi^ 10" 10"* [GcV] 

Figure 4: Ratio between the right-handed form factor of the MFV loop decay p^~'-°°P and the form factor 
of the FV tree level decay FS^ (red/full) and ratio between the MFV loop-decay form factor and the 
approximate form factor F^^ (grccn/dashcd) as function of the MFV scale /iMFV- 

with the ii — Xi mass difference ranging between 

m^^ - m^o = 9 ... 15 GeV . (54) 

We emphasize that the following results are purely illustrative. The various scenarios have not been 
required to fulfill Dark Matter constraints and/or constraints from electroweak precision data. The 
main emphasis was to achieve approximately constant masses for the NLSP and LSP. 

In Fig. H] we show, as a function of the MFV scale, the ratio of the non-resummed right-handed 
form factor i?^"^°°P in the MFV 1-loop decay to in the FV tree level decay as well as the ratio of 
to the approximate form factor F^'^^lf] As can be inferred from the figure, the approximate 
result reproduces the one-loop result down to low scales. Starting from /iMFV = 10^ GeV the finite 
terms become relevant. At /-iMFV = 10'^ GeV neglecting the finite terms in F^^^ leads to a factor 
~ 2 between the approximate and the 1-loop form factor. The non-resummed 1-loop result and the 
resummed tree level result, on the other hand, approach each other with decreasing scale of MFV 
as expected. 

Figure [5] shows the partial widths as functions of /iMFV for the approximate MFV decay, for 
the full MFV 1-loop decay and for the tree level resummed decay. An interesting feature which 
can be inferred from Fig. [5] is the size of the decay width. It does not only depend on the size of 
the logarithm but also on the coefficient of the logarithmic term, which is given in terms of the 
soft SUSY breaking parameters, particle masses and mixing angles, cf. Eq. ()43p . As explained 
above, for each value of the scale fiuFY we have chosen a different set of boundary conditions 
Mo, mi/2,Ao, tan (3, sign(//) such that the ti and Xi masses remain approximately unchanged. This 
leads for each /xmfv to a different coefficient of the logarithmic term. For hmfv = 10^^ GeV e.g. 
the parameter set and resulting masses and mixing angles are such that the coefficient becomes 

^Note that the line connecting the different points uniquely serves to guide the eye. 
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[GeV] 




IQ-l'' I 1 I 1 1 1 1 1 

10^ 10-* IC' 10** lOio 1012 IQi* lOi'' [GeV] 

^iMFV 

Figure 5: Partial decay width r{ti — >■ cxi) calculated assuming IVIFV, r^~'°°P (red/full), calculated with the 
approximate formula, F^/^ (blue/dotted), and calculated at tree level including FV, (green/dashed), 
as function of hmfv- 

rather small, so that the partial width is less than 10^^^ GeV. Due to the large value of //mfv 
the logarithmic contribution still dominates over the finite terms, however, so that there is good 
agreement between the 1-loop and approximate result. For small values of /iMFV the partial width 
can be as large as a few 10^^^ GeV as the factor, which multiplies the logarithm, turns out to be 
large for the chosen parameter set. The value of the coefficient is also the reason for the kink in 
Fig. a at fiMFY = 10^ GeV. 

Figure [5] shows, that in accordance with the behaviour of the right-handed form factors, at high 
scales the 1-loop and the approximate result agree up to the effect of the non-logarithmic terms on 
the partial width, which is at the 10% level. The 1-loop and the resummed tree level decay agree 
at low scales where the resummation effects of the large logarithms can be neglected, whereas the 
deviations are large for high scales. In summary, in order to get correct predictions for the flavour 
changing light stop decay for large scales of MFV, resummation effects have to be included. To 
further improve on this decay, the next step is the calculation of the one-loop corrections to the 
tree level stop decay including the squark mixing matrix elements from RGE evolution. This is 
deferred to a future publication. 

6 Summary and Conclusions 

In summary, we have calculated the flavour violating decay ii — )• ex? in the framework of Minimal 
Flavour Violation including also finite terms, which do not depend on the logarithm of the MFV 
scale //MFV- The one-loop decay has been compared to the approximate result derived earlier by 
Hikasa and Kobayashi which neglects the subleading terms compared to the large logarithm of 
/iMFV- It has been found that it approximates the complete one-loop result within 10% for large 
MFV scales. The approximation becomes worse with decreasing scales. The one-loop result, and 
also the approximate formula, however, do not resum the large logarithms. The resummation is 
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done by solving the renormalization group equations. Since MFV is not RG-invariant, flavour 
changing off-diagonal elements are induced in the squark mixing matrices which lead to FCNC 
couplings at tree level. They can be compared to the effective one-loop coupling in the MFV 
approach. The resummation effects turn out to be important, so that the one-loop result and 
the formula by Hikasa and Kobayashi only give an approximate value of the phenomenologically 
important light stop decay width into charm and neutralino. The next important step to improve 
the prediction for the light stop decay width will be the calculation of the one-loop corrections to 
the flavour-violating tree level decay. 



Appendix 
A Couplings 

To set up our notation for the couplings, we briefly repeat the chargino and neutralino systems. 
The chargino mass matrix, in terms of the wino mass parameter M2, fi and tan/3, is given by [30] 



where M\y denotes the charged W boson mass and we use S/3 = sin/3, = cos/3. It is diagonalized 
by two real matrices U and V, 

U*McV^' ^ U = 0. and V = { if detA^c > 

I asO^ if detMc < ^ ^ 

with the Pauli matrix a-s rendering the chargino masses positive. 0± are rotation matrices with 
the mixing angles 



tan2g_ = ^^f'"f^-^V--^) ^ 2V2MMM2S, + ,c,) _ 



The two chargino masses read 



' -\{mI + I2^ + 2M^ t [{Ml -i?f + AM^{M^c% + M| + + 2M2^iS2p)] 5 } . (58) 

ijO\ 



m ± 

x5 2 



The four-dimensional neutralino mass matrix in the {—iB, — iVFs, H^^ H2) basis has the form 



M, 



N 



I Ml -Mzswcp Mzswsp \ 

M2 Mzcwc/3 -Mzcws/3 

-MzswCj3 Mzcwcp -(U 

\ Mzswsp -Mzcwsp / 



(59) 



with c^r = 1 — s'^, = M^/M^. It can be diagonalized analytically [3T] with a single matrix Z. 

In the following, we list the couplings in the framework of MFV |3Up32I[5^ , which are needed for 
our calculation. All couplings are normalized to the weak gauge coupling g if not stated otherwise. 
Note that all charged couplings involving quarks and/or squarks have to be multiplied with the 
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CKM matrix element Vud, which we have factored out from our definition of the couphngs. 
• The couphngs of charginos and neutrahnos to the charged gauge bosons W^: 

G^^w.-G^,w^ with ^^^^^ ^ ^[Z.3C/,2 + V2Z.,t/,,] 



The couphngs of charginos and neutrahnos to charged Higgs and Goldstone bosona^i. 

M =qL,R '^fjH+ = cp[Zi^Vji + ^{Zi2 + tanewZa)Vj2] 

x'ixtH+ m+ wi ^g^^ ^ sp[Zi^U,i- j^iZ,2 + tan ewZa)Uj2] ^ ' 



r<L,R _nL,R G^G+ - S,3[Zi4Vji + ^{Zi2+tanewZii)Vj2] 

G^1x+G+ -Gi3G+ with _ I tv.^ . .^^a...r7. ^TT.^^ y°^> 



ijG+ - -cp[Zi^Uji - -^{Zi2 + taiidwZii)Uj2] 



The couphngs between neutrahnos, quarks and squarks, qi^2 — Q — x'j' 



'i2 

1 b]2 ] - V2M^ \ -so, j \ J ' ^''^ 

with ru = Zji^l sin /3 and = ^js/ cos /3 for up- and down-type quarks, and 

e% = ^/2[Z,,tw{Q,-Il) + Z,2Il] 

= -V2QgtwZji , (64) 

where tw = tan^vK- 

• The couplings between charginos, quarks and squarks, qi^2 — q' — Xj : foi' uP" ^'^'^ down-type 
(s) quarks read 



ud ( - 



aj2 



so^ j V2Mwsp 1 co^ 



hf, 1 ^ m,U,2 [ co^ 
b]i j V2Mwcp \ -so^ 



(65) 



\^ 
du 

aj2 

udu 

udu 
Oj2 




^^We work in the Fcynman gauge. 



(66) 



18 



The couplings of the gauge bosons, the charged Higgs and Goldstone bosons to quarks 
are given by 



2V2 



(67) 



rrid tan /3 + m„ cot /? „ tan /3 — rriu cot /3 



2V2M^ 2V2My, ^^^^ 

g _ -rud + ruu q _ rud + 
^•^"^ ~ 2V2MW ' ""^^ ~ 2^/2Mvy ■ 

The couphngs of the gauge bosons and charged Higgs and Goldstone bosons to squarks 
are given by 



where 



= -71" r"'^ Mf 



And for the Goldstone couplings we have 



(70) 



Cw+21 = -seu^dd Cw+22 = sensed 
Gn+ud, = ^H+ij with 
Ch+22 = se^se^su + ce^ce^S22 - sg^ce^iSi2 - ce^se^S2i (71) 

Ch+12 = -Ce^Sg^Sn+Sg^C0^S22 + Cg^Cg^Si2- Sg^S6^S21 
Ch+21 = -Se^Cg^Sii + Cg^Se^S22 - Sg^Sg^Si2 + Ce^Cg^S21 , 

Miv / . tan /3 + cot /3 



w 

S22 = ^^!^(tan/3 + cot/3) (72) 
y2Mw 

S12 = J^'^ — (/i + ^rftan/3) 
V2Mw 

S2^ = ^(f^ + ^ucotf^). 



Gg+uJ, = Canj with 

Cg+11 = ce^ce^aLL + ce^sg^aLR + SB^ce^aRL 

Cg+22 = Sg^Sg^aLL - Sg^Cg^aLR- Cg^Sg^ttRL (73) 

Cg+12 = -Cg^Sg^aiL + Cg^Cg^aiR- Sg^Sg^aRL 

Cg+21 = —Sg^Cg^QLL — Sg^Sg^aiR + Cg^Cg^aRL , 
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and 



(IRL 



w 



V2Mw 



-(/xtan/3 - 
-(/icot/3 - ^n) . 



(74) 



• For our calculation we also need the 4-squark coupling between stop, scharm and two identical 
down-type squarks. With the generation index /c = 1,2,3 and the index I = 1,2 denoting the two 
down-type squark eigenstates, it reads in terms of 



a. 



with 



mtrric 



ijdkidki 



m 



dk 



-i 



Ik 



(75) 



and 



ijll 


Vijii 






1111 








2211 






Uf. 


1211 








2111 









(76) 



The couplings with two dfc2 squarks are obtained from those with the d^i squarks by interchanging 

• Finally, the couplings between stop, scharm and two charged Higgs bosons or two charged 
Goldstone bosons in terms of are given by 

ml tan 



Gr 



2M^ 



-Q 



ijH+H+ 



(77) 



and 



where 



ml 



-Qi 



Q1IH+H+ — Ce^Cg^ Q22H+H+ — Sg^Sg^ 

Q12H+H+ = -CdtSdc Q21H+H+ = -setCdc 



(78) 



(79) 



Bl Squark self-energy contributions 



In this Appendix we give the result for the self-energy Sg^^^ in terms of the couplings defined in 
Appendix A. The self-energy receives contributions from the various diagrams in Figl2l i.e. 



^X+d ^H+d + ^G+d j.W^d 5.4,- _ 

CLtl <^Ltl C^ti CLt\ CLtl c^tl 



(80) 



20 



Note that the self-energy for ti and cr external legs is zero for vanishing c quark mass, 

t,^i^U=o = 0. (81) 

We have for 

i=l,2fc=l,..,3 ^ 

+ {af^" hf^ + af^ ^fi)^xt ^'^k ' "^x+ ' "^^^ ) } ' 

with a'^ etc. given in Eq. (j65p . The sum is taken over the chargino eigenstates and the three 
quark generations. The scalar one- loop one- and two-point integrals Ao(m) and Bo{p'^ ,mi,m2) are 
defined as [35] 



n / 2 \ ■ -4—v f 1 



(83) 
(84) 



Note the suppression by the CKM matrix elements Vcd^., Vj*^^,- We find for 5^?^^'^: 

A;=l,..,3 i=l,2 

The sum is to be taken over the three squark generations k and the two squark mass eigenstates i. 
The Goldstone contribution reads 

fc=l,..,3 1=1,2 
and the self-energy involving the W boson 

= a" E ^^'^.K.,G^.f,4.^H^+..4.{-2^o(^^) + ^o(-L) 

fc=l,..,3 i=l,2 

- (2m| + 2m2-^^ - M^)So(m| , M^., mj^^)} . (87) 
Finally, we have the tadpole contributions from the charged Higgs and Goldstone boson loop, 

^Sr(H) = 5'(-l)Ve6V;jG,-^,^^+^+Ao(m^+) + G,-^,^G+G+^o(M2,)] , (88) 
and the one from the 4-squark vertex 

= E E ^^^.^k%..4.4.^o(m^-J . (89) 

fc=l,..,3 i=l,2 
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B2 Quark self-energy 

According to the structure of the quark self-energy given in Eq. (|22p we find for the left-chiral 
contribution 

2 

Sfc(p') = g^{v,tV:,^^[-t}Bi{p',mt,mH+)-B,{p',mb,Mw)] 

i=l,2A;=l,..,3 j=l,2 

+ E ^V,,,V;,J-l-2i?i(p2,md„M^)]}, (90) 

A;=l,..,3 

with cLji'^iCL^'l^ given in Eq. (|66p . The sums are taken over all possible chargino states (j = 1,2), 
the three quark and squark generations {k = 1, 2, 3) and the two squark mass eigenstates (i = 1,2). 
Furthermore, Bi in terms of the scalar one- and two-point functions is given by 

Si(p^,mo,mi) = ^ [^o("T-o) -^o("ii) - (P^ -"i? + m^)So(p^,mo,mi)] . (91) 

The right-chiral contribution reads 

Sfe(p') = E Ved,V,\^[-cot2Si(/,mrf,,mH+)-Si(p2,m,,,MH^)] 

+ E E E^^'^.^*<^.(-^i^&^f)5l(p^-,;,-J^^ (92) 

i=l,2fc=l,2,3 4=1,2 

with h'^'l'^ ,})^^^ defined in Eq. ()66p . It vanishes for rric = 0. The left-chiral scalar contribution can 
be cast into the form 

2 

+ E E E^^.^.-x;4'^t^^o(/, -,;)}, (93) 

j=l,2fc=l,..,3i=l,2 

which also vanishes for zero charm quark mass. For the right-chiral scalar contribution we find 



^fcip") = gHv,,V:J^[Bo{p\mH^,m,)-B^{p\Mw,m,)] 

+ E E E^^^.K..-,.6J^*a^^i?o(p^-4.m,.)}. (94) 

i=l,2fc=l,..,3 j=l,2 

Note that in case of real two-point functions we have Ti^^ = S^/*. 
B3 Vertex correction 

The vertex contributions to the left-chiral form factor vanish for rric = 0. For the right-chiral 
form factor they are given by the various right-chiral contributions from the vertex correction 
graphs depicted in Figl2] 

= i [r^+Jd + ^x+H+d + ^x+G+d + ^x+W+d + ^dH+d + ^dG+d + ^x+dH+ + ^x+dG+ 

+'^MW^+'^dx+wAi^l)- (95) 
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We have for 

-C2ijfc(m^+m|o - m?^mj^p + {-c^ijkm^omd^ + C4ijkm-+md^){m^^ - m|o)] 

Co (m| , m|o , 0, m -+ , m^^ , m^^^ ) | / (m| - m|o ) , (96) 

where the sum over all possible chargino eigenstates (j = 1,2), all three generations of down 
type quarks and squarks (A; = 1,2,3) as well as the two squark mass eigenstates {i = 1,2) has to 
be taken. We have introduced the abbreviations 

di.c tdi^ dh. I dh-c-itdi^ du di^citdi^-idu 

Clijk = aj'^ aj^''a^'irn^++aji b-^^a^^ma^ C2ijk = a-1 h-^^h^\ ^^^^ 

duCitdu dh dhC tdu idh 

C3ijk = flj- ai- C4ijk = a^{h^^ , 

with the various couplings defined in Eqs. (j63l65l66p . The scalar one-loop 3-point function is given 
by 

Co{pI,pI,{Pi +P2f, mi, 1712,1713,) = 

_^ ^ f d Jxi 1 , , 

"'^ J (2^)" (fc2 - mf) [(fc + pi)2 - ml] [{k + + p2)2 - mj] ' ^^^^ 

We find 

^x+H+dim'^iJ = / Yl ^cbVtb ^^^^^^ {[cijm^o+C2jm'ljBo{m'l^,m^+,mb) 

j=l,2 \J IMy^ J 

- [cij-m^o + C2jm|o]-Bo(m|o, mj|^+, m-+) + [cijm^o(m| - m|o - m-6 + m^+) 
+C2j(m|m|-+ - m^^oml) + m^+{c3jmb - C4jm^o)(m| - m|o)] 

Co (m| , m|o , 0, mfo, m -+ , mH+ ) } / (m| - m|o ) , (99) 

with 

cij = a%G^.^+mh + b%Gfj4m^+ C2j = b%Gf^^+ 

csj = afiG^.^+ C4j = hfiG[.^+ . 

Note that we set the down and strange quark mass to zero, md = mg = 0. We have 

VG+.(m|) = V,,K,^^{[c5,m,o + C6,m?jSo(m|,m-.,m,) 

j=l 2 V 2,IVl\Y J 

- [c5j?7i-o + C6jm|o]5o(m|o , M^/, m-+ ) + [cs^m^o (m| - m|o - "i^ + M^) 
+C6j(m|M^ - mioml) + m^+{cjjmi, - C8jm^o)(m| - m|o)] 

Co(m|,m|o,0,mb,m-+,MH/)|/(m| - m|o) , (100) 



23 



where 

C7j = af^G^.(j+ csj = bf^G^.Q+ . 

And 

j=l,2fc=1...3 ' 
+ 2[cijfcm^o + C2jk{m^^ - m|o)]-Bo(m|o, M^y, m-+) + 2c2jfc(m| - m|o) 

Bo{0, Mw, maj + [2c2jk{inl^ + - ?n|)(m| - m|o) + 2m^ocijk 
(m|o - m| + m^^ - M^) + 2m-+(c3jfcm^o + 2cijkmdJ{mf^ - m|o)] 
Co(m|,m|o,0,md^,m-+,M,4/) - 2c2jfc(m| -m|o)|/(m?^ ""^1°) ' (^0^) 



with 



Next 



and 



[a\,m,{ml - m|o) - b^.m^oK - m^+)]Co(m|, m|o, 0, , m^^, m,)} 
/(m|-m|o) (103) 



i=i,2 v2Aivy 
[aijmfe(m| -m|o) - h\^m^o{ml - M^)]Co(m|, m|o, 0, M^y, m^^, mf,)| 



Furthermore, 



j=l,2fc=1...3i=l,2 

Co(m|,m|o,0,mj^^,m^+,m-+)|/(m| - m?,,) (105) 



with 



c^ijk = GvjH+^i" ^Sijfc = G^-H+afi" , (106) 
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and 



with 



We have 



^ddW+ 



(m? 



And finally 



+ 



with 



Co (m| , m|o , , m , Mw , ) } / (' 



^X+dG+("^h) ~ f ^ X] X] X] ^G+hdk 
jr=l,2fe=1...3i=l,2 



+ Bo{m'lo,m-+,Mw)] + 



iml-mio) , 



C7ijk 



= G 



■R 



df^c 



(107) 



(108) 



= E , E ^ E„ ^^^Ve., Vr,, { [-ail (ml + m|o) - ^?>,om,J 



j=l,2fc=1...3't=l,2 

Bo (m| , m , , Miy ) + [2a'^J^ m|o + ^^^^ m^o ma, ] Bo (m|o , , m , ) 



+ i^l - mlo)Bo{0, ma, , Mw) + [2^^ (m|o (m|o - m| - - ^mj^ + M^) 

+m|(mj^^ - ^mjj) - 6^*^m^omrf^ (2m| + mj^ - M 
Co (m| , m|o , 0, Mh/ , m^^^ , nid^ ) | / (m| - m|o ) • 



■ 2m|o)] 



(109) 



■^d^VcdXtdk { hiifc^l + cioijfcm^om-+]5o("i| , Mw,m^J 



E E E 

j=l,2A:=1...3j=l,2 

[c9ijfc(m'|o - 2m? ) - cioijfcm^om-+]Bo(m|o, Mvi/) + [cQijkimj^m^-.-^ 

, 9.9 9.9m 



^1 

Co (m| , m|o , 0, m , Miy , m -+ ) | / (m| - m|o ) , 



2 _2 



GL df^c 



— 2mr m T 
(110) 

(111) 



C FCNC counterterm 



We start from the u — u— neutralino part of the Lagrangian in the interaction basis, expressed in 
terms of the bare squark and quark fields, uf'^ and uf^^ and the bare quark mass matrix T^fj\ 
where i,j = 1,2,3 denote the generation indices and I = 1, ..,4 the neutralino mass eigenstates, 



-(0) 



r,0 



,(0) 



9 e-Ri UiR '^L Xi + 



gZumf \ (0) 
\/2Mv.sin/3j"^-^^^^' 

V2Mwsm/3 



uJlVLX'i + h.c. 



(112) 
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The couplings e^jn have been defined in Eq. (I64p . Let us look at the right-chiral part of the 
coupling. Rotation to the mass eigenstates yields 



mux" 



-n7 rr.(°) 

m{0)jjUL{0) ^Ui TT/(0)t ~m(0) -p ^0 , -m(0) ..n^ (0) ^ff^m^ (o)t ^(0) t, ^0 , 



\/2MH/sin/3 



{i,j,k = 1,2,3, s = 1,..,6) 
Note, that W^f^^ = wf^^, wf^l^ = cf. Eq. Upon renormahzation we replace [M 



c. 

(113) 



+ 



2 -J 



^{0)t~™(0) 



(1 + 5n"^)[7"^ 



m 



(0) 



m + 5m , 



(114) 

(115) 
(116) 



where we have suppressed the indices. The wave function denotes a six-component column 
vector. With the replacement Wl,r = Wl,rU^^'^ , cf. Eq. (fT5]) . we have for the Yukawa part of 
the coupling 

1 + j (1 + (^'u"^)[/"^(m + (5m)C/"«t(i + 5u"«t)VFjj(l + 5u)t) h + _ j , (117) 

times {—gZi/i)/[^j2Mwsp). For the mass renormahzation we choose the renormahzation prescrip- 
tion such that the bare mass matrices and hence 6m are diagonal, i.e. 



(1 + (5n"^)C/"^(m + (^m)C/"«t(i + j^^^flt) = + 



(118) 



where I? denotes diagonal matrices. This is possible since the off-diagonal elements can be ab- 
sorbed into the off-diagonal elements of the antihermitian part of the right-handed wave function 
renormahzation matrices [36] ■ Exploiting the unitary of the mixing matrices we finally find for the 
renormalized Lagrangian in the mass eigenstate basis 



with the couplings given by 



^isl — 




^isl — 




SG^si = 




2 




gZiA 



fw sin /3 



(119) 

(120) 
(121) 



\/2MvKsin/3 



(122) 
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5G^ 



isl 



Uj 



9Z14 



ts 



\f2Mw sin P 



+muAj{wl)jt6wlg + 5muAj(wl)js 



(123) 



In the framework of MFV at /Umfv the W matrix is diagonal in flavour space at tree level. At one- 
loop flavour off-diagonal elements are induced through the wave function and the mixing matrix 
renormalization. 
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